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Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically in- 
volve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, 
among other things, on the separation between the two detectors: the smaller the separation, the better the 
sensitivity. Hence, a co-located detector pair is more sensitive to a gravitational-wave background than a non- 
co-located detector pair. However, co-located detectors are also expected to suffer from correlated noise from 
instrumental and environmental effects that could contaminate the measurement of the background. Hence, 
methods to identify and mitigate the effects of correlated noise are necessary to achieve the potential increase 
in sensitivity of co-located detectors. Here we report on the first SGWB analysis using the two LIGO Hanford 
detectors and address the complications arising from correlated environmental noise. We apply correlated noise 
identification and mitigation techniques to data taken by the two LIGO Hanford detectors, HI and H2, during 
LIGO’s fifth science run. At low frequencies, 40 — 460 Hz, we are unable to sufhciently mitigate the correlated 
noise to a level where we may confidently measure or bound the stochastic gravitational-wave signal. However, 
at high frequencies, 460 — 1000 Hz, these techniques are sufficient to set a 95% confidence level (C.L.) upper 
limit on the gravitational-wave energy density of f2(/) < 7.7 x 10~ 4 (//900 Hz) 3 , which improves on the 
previous upper limit by a factor of ~ 180. In doing so, we demonstrate techniques that will be useful for future 
searches using advanced detectors, where correlated noise (e.g., from global magnetic fields) may affect even 
widely separated detectors. 


I. INTRODUCTION 

The detection of a stochastic gravitational-wave back- 
ground (SGWB), of either cosmological or astrophysical ori- 
gin, is a major science goal for both current and planned 
searches for gravitational waves (GWs) M- Given the 
weakness of the gravitational interaction, cosmological GWs 
are expected to decouple from matter in the early universe 
much earlier than any other form of radiation (e.g., photons, 
neutrinos, etc.). The detection of such a primordial GW back- 
ground by the current ground-based detectors , proposed 
space-based detectors |8|,|9|], or a pulsar timing array lUOlnin 
would give us a picture of the universe mere fractions of a 
second after the Big-Bang GHSE3, allowing us to study the 
physics of the highest energy scales, unachievable in standard 
laboratory experiments 01 . The recent results from the BI- 
CEP2 experiment indicate the existence of cosmic microwave 
background B-mode polarization at degree angular scales JH, 
which may be due to an ultra-low frequency primordial GW 
background, such as would be generated by amplification of 
vacuum fluctuations during cosmological inflation; however, 
it cannot currently be ruled out that the observed B-mode po- 
larization is due to a Galactic dust foreground llf4[l~i~5ll ). These 
GWs and their high frequency counterparts in standard slow- 
roll inflationary model are several orders of magnitude below 
the sensitivity levels of current and advanced LIGO detectors. 
Hence they are not the target of our current analysis. How- 
ever, many non-standard inflationary models predict GWs that 
could be detected by advanced LIGO detectors. 

On the other hand, the detection of a SGWB due to spatially 
and temporally unresolved foreground astrophysical sources 
such as magnetars 0, rotating neutron stars 0, galactic 
and extragalactic compact binaries lfl8l420ll . or the inspiral and 


collisions of supermassive black holes associated with distant 
galaxy mergers m ], would provide information about the spa- 
tial distribution and formation rate of these various source pop- 
ulations. 

Given the random nature of a SGWB, searches require 
cross-correlating data from two or more detectors 1 11 122142511 . 
under the assumption that correlated noise between any two 
detectors is negligible. For such a case, the contribution to 
the cross-correlation from the (common) GW signal grows lin- 
early with the observation time T, while that from the noise 
grows like s/T. Thus, the signal-to-noise ratio (SNR) also 
grows like s/T. This allows one to search for stochastic sig- 
nals buried within the detector noise by integrating for a suffi- 
ciently long interval of time. 

For the widely-separated detectors in Livingston, LA and 
Hanford, WA, the physical separation (~ 3000 km) elimi- 
nates the coupling of local instrumental and environmental 
noise between the two detectors, while global disturbances 
such as electromagnetic resonances are at a sufficiently low 
level that they are not observable in coherence measurements 
between the ( first-generation) detectors at their design sensi- 
tivity immu. 

While physically-separated detectors have the advantage of 
reduced correlated noise, they have the disadvantage of re- 
duced sensitivity to a SGWB; physically-separated detectors 
respond at different times to GWs from different directions 
and with differing response amplitudes depending on the rel- 
ative orientation and (mis)alignment of the detectors li23l - i25ll . 
Co-located and co-aligned detectors, on the other hand, such 
as the 4 km and 2 km interferometers in Hanford, WA (de- 
noted HI and H2), respond identically to GWs from all di- 
rections and for all frequencies below a few kHz. They are 
thus, potentially, an order-of-magnitude more sensitive to a 
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SGWB than e.g., the Hanford-Livingston LIGO pair. But this 
potential gain in sensitivity can be offset by the presence of 
correlated instrumental and environmental noise, given that 
the two detectors share the same local environment. Meth- 
ods to identify and mitigate the effects of correlated noise are 
thus needed to realize the potential increase in sensitivity of 
co-located detectors. 

In this paper, we apply several noise identification and mit- 
igation techniques to data taken by the two LIGO Hanford 
detectors, HI and H2, during LIGO’s fifth science run (S5, 
November 4, 2005, to September 30, 2007) in the context of 
a search for a SGWB. This is the first stochastic analysis us- 
ing LIGO science data that addresses the complications intro- 
duced by correlated environmental noise. As discussed in the 
references lf29i[30il , the coupling of global magnetic fields to 
non-colocated advanced LIGO detectors could produce signif- 
icant correlations between them thereby reducing their sensi- 
tivity to SGWB by an order of magnitude. We expect the cur- 
rent H1-H2 analysis to provide a useful precedent for SGWB 
searches with advanced detectors in such (expected) corre- 
lated noise environment. 

Results are presented at different stages of cleaning applied 
to the data. We split the analysis into two parts — one for 
the frequency band 460-1000 Hz, where we are able to suc- 
cessfully identify and exclude significant narrow-band correla- 
tions; and the other for the band 80-160 Hz, where even after 
applying the noise reduction methods there is still evidence 
of residual contamination, resulting in a large systematic un- 
certainty for this band. The frequencies below 80 Hz and be- 
tween 160-460 Hz are not included in the analysis because of 
poor detector sensitivity and contamination by known noise 
artifacts. We observe no evidence of a SGWB and so our final 
results are given in the form of upper-limits. Due to the pres- 
ence of residual correlated noise between 80-160 Hz, we do 
not set any upper-limit for this frequency band. Since we do 
not observe any such residual noise between 460-1000 Hz, in 
that frequency band and the 5 sub-bands assigned to it, we set 
astrophysical upper-limits on the energy density of stochastic 
GWs. 

The rest of the paper is organized as follows. In Sec. [Til we 
describe sources of correlated noise in HI and H2, and the en- 
vironmental and instrumental monitoring system. In Sec. [Ill] 
we describe the cross-correlation procedure used to search for 
a SGWB. In Secs.|IV]and[V]we describe the methods that we 
used to identify correlated noise, and the steps that we took 
to mitigate it. In Secs. I VII and I VIII we give the results of our 
analysis applied to the S5 H1-H2 data. Finally, in Sec. lVIIIl 
we summarize our results and discuss potential improvements 
to the methods discussed in this paper. 


II. COMMON NOISE IN THE TWO LIGO HANFORD 
DETECTORS 

At each of the LIGO observatory sites the detectors are 
supplemented with a set of sensors to monitor the local envi- 
ronment SEl. Seismometers and accelerometers measure 
vibrations of the ground and various detector components; 


microphones monitor acoustic noise; magnetometers monitor 
magnetic fields that could couple to the test masses (end mir- 
rors of the interferometers) via the magnets attached to the 
test masses to control their positions; radio receivers monitor 
radio frequency (RF) power around the laser modulation fre- 
quencies, and voltage line monitors record fluctuations in the 
AC power. These physical environment monitoring (PEM) 
channels are used to detect instrumental and environmental 
disturbances that can couple to the GW strain channel. We 
assume that these channels are completely insensitive to GW 
strain. The PEM channels are placed at strategic locations 
around the observatory, especially near the corner and ends 
of the L-shaped interferometer where important laser, optical, 
and suspension systems reside in addition to the test masses 
themselves. 

Information provided by the PEM channels is used in many 
different ways. The most basic application is the creation 
of numerous data quality flags identifying stretches of data 
that are corrupted by instrumental or environmental noise 0. 
The signals from PEM channels are critical in defining these 
flags; microphones register airplanes flying overhead, seis- 
mometers and accelerometers detect elevated seismic activity 
or anthropogenic events (trucks, trains, logging), and magne- 
tometers detect fluctuations in the mains power supply and the 
Earth’s magnetic field. 

In searches for transient GW signals, such as burst or coa- 
lescing binary events, information from the PEM channels has 
been used to construct vetoes J3314361] . When a clear associ- 
ation can be made between a measured environmental event 
and a coincident glitch in the output channel of the detector, 
then these times are excluded from the transient GW searches. 
These event-by-event vetoes exclude times of order hundreds 
of milliseconds to a few seconds. 

Similarly, noise at specific frequencies, called noise lines, 
can affect searches for GWs from rotating neutron stars or 
even for a SGWB. In S5, data from PEM channels were used 
to verify that some of the apparent periodic signals were in 
fact due to noise sources at the observatories [37- 3*8 1. Typ- 
ically the neutron-star search algorithms can also be applied 
to the PEM data to find channels that have noise lines at the 
same frequencies as those in the detector output channel. The 
coherence is also calculated between the detector output and 
the PEM channels, and these results provide additional infor- 
mation for determining the source of noise lines. 

The study of noise lines has also benefited past LIGO 
searches for stochastic GWs. For example, in LIGO’s search 
for a SGWB using the data from the S4 run l27ll . correlated 
noise between the Hanford and Livingston detectors was ob- 
served in the form of a forest of sharp 1 Hz harmonic lines. It 
was subsequently determined that these lines were caused by 
the sharp ramp of a one-pulse-per-second signal, injected into 
the data acquisition system to synchronize it with the Global 
Positioning System (GPS) time reference. In the S5 stochastic 
search (280, there were other prominent noise lines that were 
subsequently identified through the use of the PEM signals. 

In addition to passive studies, where the PEM signals are 
observed and associations are made to detector noise, there 
have also been a series of active investigations where noise 
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was injected into the detector environment in order to measure 
its coupling to the GW channel. Acoustic, seismic, magnetic, 
and RF electromagnetic noise were injected into the observa- 
tory environment at various locations and responses of the de- 
tectors were studied. These tests provided clues and ways to 
better isolate the detectors from the environment. 

All the previous LIGO searches for a SGWB have used the 
physically-separated Hanford and Livingston detectors and as- 
sumed that common noise between these non-colocated de- 
tectors was inconsequential. This assumption was strongly 
supported by observations — i.e., none of the coherence mea- 
surements performed to date between these detectors revealed 
the presence of correlations other than those known to be in- 
troduced by the instrument itself (for example, harmonics of 
the 60 Hz power line). Since the analysis presented here uses 
the two co-located Hanford detectors, which are susceptible 
to correlated noise due to the local environment, new methods 
were required to identify and mitigate the correlated noise. 


III. CROSS-CORRELATION PROCEDURE 


The energy density spectrum of SGWB is defined as 

(1) 

Pc df 

3 c 2 H 2 

where p c (= 87rG ° ) is the critical energy density and p gw is 
the GW energy density contained in the frequency range / 
and / + df. Since most theoretical models of stochastic back- 
grounds in the LIGO band are characterized by a power-law 
spectrum, we will assume that the fractional energy density in 
GWs |39l| has the form 

figw(/) = 0, , (2) 

where a is the spectral index and / re f is some reference fre- 
quency. We will consider two values for the spectral index: 
a = 0 which is representative of many cosmological mod- 
els, and a = 3 which is characteristic of many astrophysical 
models. This latter case corresponds to a flat (i.e., constant) 
one-sided power spectral density (PSD) in the strain output of 
a detector 5 gw (/), since 


Sgw(/) 


3 Hi O gw (/) 

107T 2 f 3 


( 3 ) 


Here Hq is the present value of the Hubble parameter, as- 
sumed to be Hq = 68 km/s/Mpc 

Following the procedures described in |25], we construct 
our cross-correlation statistic as estimators of fl Q for individ- 
ual frequency bins, of width A/, centered at each (positive) 
frequency /. These estimators are simply the measured val- 
ues of the cross-spectrum of the strain output of two detectors 
divided by the expected shape of the cross-correlation due to 
a GW background with spectral index a: 




2 TO(/)] 

T 7 (f)S a (f) 


( 4 ) 


Here T is the duration of the data segments used for Fourier 
transforms; Si(/), $ 2 (f) are the Fourier transforms of the 
strain time-series in the two detectors; S a (f) is proportional 
to the assumed spectral shape. 


s a (f) = 


3 H 2 q 1 

107T 2 f 3 



( 5 ) 


and 7 (/) is the overlap reduction function l23l - l25ll . which en- 
codes the reduction in sensitivity due to the separation and rel- 
ative alignment of the two detectors. For the H1-H2 detector 
pair, 7 (/) 1 for all frequencies below a few kHz(41j. 

In the absence of correlated noise, one can show that 
the above estimators are optimal — i.e., they are unbiased, 
minimal-variance estimators of S1 Q for stochastic background 
signals with spectral index a. Assuming that the detector 
noise is Gaussian, stationary, and much larger in magnitude 
than the GW signal, the expectation value of the variance of 
the estimators is given by 


1 Pi(f)P 2 (f) 
2TA/7 2 (/)S 2 (/)’ 


(6) 


where Pi(f), Pi{f) are the one-sided PSDs of the detector 
output si (/), §2 (/) respectively. For a frequency band con- 
sisting of several bins of width A/, the optimal estimator and 
corresponding variance are given by the weighted sum 
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( 7 ) 


A similar weighted sum can be used to optimally combine the 
estimators calculated for different time intervals [42]. 

In the presence of correlated noise, the estimators are bi- 
ased. The expected values are then 

(&a{f)) = V> a +ri a {f), ( 8 ) 


where 


vM) = 


®[Nl2(f)] 

7 (f)Sa(f) ■ 


( 9 ) 


Here A/^12 (/) = jt(n*(f)n 2 {f)) is the one-sided cross- 
spectral density (CSD) of the correlated noise contribution 
to Si, s 2 - The expression for the variance er 2 (/) 
is unchanged in the presence of correlated noise provided 
| iVi2 (/) | ■C Pi(f),p 2 (f). For the summed estimator Cl a , 
we have 

(Ci a ) = n a + p a (io) 


where 


Va = 




(ii) 


is the contribution from correlated noise averaged over time 
(not shown) and frequency. Thus, correlated noise biases our 
estimates of the amplitude of a SGWB. Here we also note that 
rj a can be positive or negative while (>,, is positive by defi- 
nition. The purpose of the noise identification and removal 
methods that we describe below is to reduce this bias as much 
as possible. 
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IV. METHODS FOR IDENTIFYING CORRELATED NOISE 
A. Coherence calculation 


Perhaps the simplest method for identifying correlated 
noise in the H1-H2 data is to calculate the magnitude squared 
coherence, ri 2 (/) = |7i2(/)| 2 , where 


_ 2 (5*(/)5 2 (/))iy 


( 12 ) 


Here T denotes the duration of a single segment of data, and 
angle brackets ( ) jv denotes an average over N segments used 
to estimate the CSD and PSDs that enter the expression for 
7i2- If there are no correlations (either due to noise or a GW 
a signal) in the data, the expected value of T 12 (f) is equal to 
1 /N. This method is especially useful at finding narrowband 
features that stick out above the expected 1 /N level. Since we 
expect a SGWB to be broadband, with relatively little varia- 
tion in the LIGO band (^80-1000 Hz), most of these features 
can be attributed to instrumental and/or environmental corre- 
lations. We further investigate these lines with data from other 
PEM channels and once we confirm that they are indeed en- 
vironmental/instrumental artifacts, we remove them from our 
analysis. 

Plots of Ti2 (/) for three different frequency resolutions 
are shown in Figs. |T] and [2] for two frequency bands, 80- 
160 Hz and 460-1000 Hz, respectively. In Fig. Q] note the 
relatively wide structure around 120 Hz, which is especially 
prominent in the bottom panel where the frequency resolu- 
tion is 100 mHz. This structure arises from low-frequency 
noise (dominated by seismic and other mechanical noise) up- 
converting to frequencies around the 60 Hz harmonics via a 
bilinear coupling mechanism. While these coupling mecha- 
nisms are not fully understood, we reject the band from 102- 
126 Hz for our analysis, given the elevated correlated noise 
seen in this band. (A similar plot at slightly lower and higher 
frequencies shows similar noisy bands from 40-80 Hz and 
160-200 Hz.) A closer look at the coherence also identifies 
smaller structures at 86-90 Hz, 100 Hz, 140-141 Hz, and 
150 Hz. A follow-up analysis of PEM channels (which is dis- 
cussed in more detail later) revealed that the grayed bands in 
Figs. |T|and [2] were highly contaminated with acoustic noise or 
by low-frequency seismic noise up-converting to frequencies 
around the 60 Hz harmonics via a bilinear coupling mecha- 
nism; so we rejected these frequency bands from subsequent 
analysis. As mentioned earlier, the 1 60-460 Hz band was 
not used in this analysis, because of similar acoustic and seis- 
mic contamination, as well as violin-mode resonances of the 
mirror-suspension wires (see Sec. lIVDt . 

As shown in Fig. [2] the coherence at high frequencies (460- 
1000 Hz) is relatively clean. The only evidence of narrow- 
band correlated noise is in ±2 Hz bands around the 60-Hz 
power-line harmonics, and violin-mode resonances of mirror 
suspensions at 688.5 ± 2.8 Hz and 697 ± 3.1 Hz. The ele- 
vated coherence near 750 Hz at 100 mHz resolution is due 
to acoustic noise coupling to the GW channels. Notching the 
power-line harmonics and violin-mode resonances amounts to 
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FIG. 1. Coherence P 12 between HI and H2 computed in the fre- 
quency band 80-160 Hz using all of the S5 data, for three differ- 
ent frequency resolutions: 1 mHz, 10 mHz, and 100 mHz (top-to- 
bottom). The insets show that the histograms of the coherence at the 
analyzed frequencies follow the expected exponential distribution for 
Gaussian noise, as well as the presence of a long tail of high coher- 
ence values at notched frequencies. A stochastic broadband GW sig- 
nal of SNR = 5 would appear at a level of < 10 x below the dashed 
l/N line. 
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FIG. 2. Same as Fig. |T| but at higher frequencies, 460-1000 Hz. 
Note the coherence peaks at the harmonics of the 60 Hz power lines 
(notched in the analysis). The elevated coherence near 750 Hz at 
100 mHz resolution is due to acoustic noise coupling to the GW chan- 
nels. The long tail in the 100 mHz plot is due to excess noise around 
750 Hz, which was removed from the final analysis using PEM notch- 
ings (see Sec. 0. A stochastic broadband GW signal of SNR = 5 
would appear at a level of < 10 x below the dashed 1/JV line. 


the removal of ~ 9% of the frequency bins over the entire 
high-frequency band. 


B. Time-shift analysis 

A second method for identifying narrowband correlated 
noise is to time-shift the time-series output of one detec- 
tor relative to that of the other detector before doing the 
cross-correlation analysis (430. By introducing a shift of 
±1 second, which is significantly larger than the correla- 
tion time for a broadband GW signal (~ 10 ms, cf. Fig. [9]), 
we eliminate broadband GW correlations while preserving 
narrowband noise features. Using segments of duration 
T = 1 s, we calculate the time-shifted estimators Cl a ,TS (/)> 
variance c>q TS (/), and their ratio SNRq c ,,TS (/) = 

^a,Ts(/)/o'o a .Ts(/)- The calibration and conditioning of 
the data is performed in exactly the same way as for the final 
search, which is described in detail in Secs.IVland lVIl 

We excise any frequency bin with ISNRq^ ,ts (/) | > 2 on 
the grounds that it is likely contaminated by correlated noise. 
This threshold was chosen on the basis of initial studies per- 
formed using playground data to understand the effectiveness 
of such cut. This criterion can be checked for different time- 
scales, such as weeks, months, or the entire data set. This 
allows us to identify transient effects on different time-scales, 
which may be diluted (and unobservable) when averaged over 
the entire data set. 


C. PEM coherence calculations 


Another method for identifying correlated noise is to first 
try to identify the noise sources that couple into the individ- 
ual detector outputs by calculating the coherence of si and s 2 
with various PEM channels zf. 


= 2 (Sj(f)zi(f))N 

~ T vW/)MW)>jv ' 


(13) 


Here *=1,2 labels the detector outputs and I labels the PEM 
channels. For our analysis we used 172 PEM channels located 
near the two detectors. In addition to the PEM channels, we 
used a couple of auxiliary channels associated with the stabi- 
lization of the frequency of the lasers used in the detectors, 
which potentially carry information about instrumental corre- 
lations between the two detectors. (Hereafter, the usage of the 
acronym PEM will also include these two auxiliary channels.) 
The Fourier transforms are calculated for each minute of data 
(T = 60 s), and the average CSDs and PSDs are computed 
for extended time-periods — weeks, months, or the entire run. 
We then perform the following maximization over all PEM 
channels, for each frequency bin /, defining: 


7!2,pem(/) = maxSR [ 7 i/(/) x 7 21 (f)] • (14) 


Note that by construction 7i2,pem(/) is real. 
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As discussed in 0 , 7i2,pem(/) is an estimate of the in- 
strumental or environmental contribution to the coherence be- 
tween the GW channels of HI and H2. This estimate is only 
approximate, however, and potentially suffers from systematic 
errors for a few reasons. First, the PEM coverage of the obser- 
vatory may be incomplete — i.e., there may be environmental 
or instrumental effects that are not captured by the existing 
array of PEMs. Second, some of the PEM channels may be 
correlated. Hence, a rigorous approach would require calcu- 
lating a matrix of elements 7 jj(/), and then inverting this 
matrix or solving a set of linear equations involving elements 
°f 7 jj(/). In practice, due to the large number of channels 
and the large amount of data, this is a formidable task. In- 
stead, we simply maximize, frequency-by-frequency, over the 
contributions from different PEM channels and use this max- 
imum as an estimate of the overall environmental contribu- 
tion to 712 (/). Finally, these coherence methods do not take 
into account the nonlinear upconversion processes in which 
low-frequency disturbances, primarily seismic activity, excite 
higher-frequency modes in the instrument. 

Since the measured signal-to-noise ratio for the estimator 
can be written as 


SNR(/) = V2TA/K [7i 2 (/)] , (15) 

we can simply approximate the contribution of the PEM chan- 
nels to the stochastic GW signal-to-noise ratio as 


SNRpem(/) = \/ 2TAf 712, pem (/) , (16) 

remembering that 7i2,pem(/) is real. The PEM contribution 
to the estimators ff Q (/) is then 

fl a ,PEM(/) = SNRpem(/)cto q (/) (17) 

where a 7 (/) is the statistical uncertainty defined by Eq.[ 6 ] 

We can use the PEM coherence calculations in two comple- 
mentary ways. First, we can identify frequency bins with par- 
ticularly large instrumental or environmental contributions by 
placing a threshold on |SNRpem(/)| and exclude them from 
the analysis. Second, the frequency bins that pass this data- 
quality cut may still contain some residual environmental con- 
tamination. We can estimate at least part of this residual con- 
tamination by using f/ a ,PEM(/) for the remaining frequency 
bins. 

As part of the analysis procedure, we were able to identify 
the PEM channels that were responsible for the largest coher- 
ent noise between the GW channels in HI and H2 for each 
frequency bin. For both the low and high frequency analyses, 
microphones and accelerometers in the central building near 
the beam splitters of each interferometer registered the most 
significant noise. Within approximately 1 Hz of the 60-Hz har- 
monics, magnetometers and voltage line monitors registered 
the largest correlated noise, but these frequencies were already 
removed from the analysis due to the significant coherence 
(noise) level at these frequencies, as mentioned in Sec. lIV Al 
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FIG. 3. (Color online) Comparison of the (absolute value of the) 
SNRs calculated by the PEM-coherence and the time-shift tech- 
niques. The vertical dotted lines indicate the frequency bands used 
for the low (80-160 Hz; black dotted lines) and high (460-1000 Hz; 
magenta dotted lines) frequency analyses. Note that SNRn^ ,ts (/) 
is a true signal-to-noise ratio, so values < 2 are dominated by ran- 
dom statistical fluctuations. SNRq Qi pem(/), on the other hand, is 
an estimate of the PEM contribution to the signal-to-noise ratio, so 
values even much lower than 2 are meaningful measurements (i.e., 
they are not statistical fluctuations). The two methods agree very 
well in identifying contaminated frequency bins or bands. Note that 
both methods indicate that the 80-160 Hz and 460-1000 Hz bands 
have relatively low levels of contamination. 


D. Comparing PEM-coherence and time-shift methods 


Figure[3]shows a comparison of the SNRs calculated by the 
PEM-coherence and time-shift methods. The agreement be- 
tween these two very different techniques in identifying con- 
taminated frequency bins (those with |£>iV.R| > a few) is re- 
markably good, which is an indication of their robustness and 
effectiveness. Moreover, Fig. [3] shows that the frequency re- 
gion between 200 Hz and 460 Hz is particularly contaminated 
by environmental and/or instrumental effects. Hence, in this 
analysis we focus on the low-frequency region (80-160 Hz) 
which is the most sensitive to cosmological backgrounds (i.e., 
spectral index a = 0 ), and on the high-frequency region 
(460-1000 Hz) which is less contaminated and more suitable 
for searches for astrophysically-generated backgrounds (e.g., 
a = 3). 

We emphasize that the PEM channels only monitor the in- 
strument and the environment, and are not sensitive to GWs. 
Similarly, the time-shift analysis, with a time-shift of ±1 sec- 
ond, is insensitive to broad-band GW signals. Hence, any data- 
quality cuts based on the PEM and time-shift studies will not 
affect the astrophysical signatures in the data — i.e., they do 
not bias our estimates of the amplitude of a SGWB. 
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E. Other potential non-astrophysical sources of correlation 

We note that any correlations that are produced by envi- 
ronmental signals that are not detected by the PEM sensors 
will not be detected by the PEM-coherence technique. Fur- 
thermore, if such correlations, or correlations from a non- 
environmental source, are broadband and flat (i.e., do not vary 
with frequency over our band), they will not be detected by 
either the PEM-coherence or the time-shift method. One po- 
tential source of broadband correlation between the two GW 
channels is the data acquisition system itself. We investigated 
this possibility by looking for correlations between 153 chan- 
nel pairs that had no physical reason to be correlated. We 
found no broadband correlations, although we did find an un- 
explained narrow-band correlation at 281.5 Hz between 10 of 
153 channel pairs. Note that 281.5 Hz is outside of the fre- 
quency bands analyzed in this study. 

We addressed the potential of correlations from un- 
monitored environmental signals by searching for coupling 
sites four times over the course of the run by injecting large 
but localized acoustic, seismic, magnetic and RF signals. New 
sensors were installed at the two coupling sites that had the 
least coverage. However, we found that the new sensors, even 
after scaling up to the full analysis period, contribute less than 
1% of the total frequency notches; hence it is safe to assume 
that we had sufficient PEM coverage throughout our analysis 
period. 

We also examined the possibility of correlations between 
the HI and H2 detectors being generated by scattered light. 
We considered two mechanisms: first, light scattered from one 
detector affecting the other detector, and second, light from 
both detectors scattering off of the same site and returning 
to the originating detectors. We did not observe, and do not 
expect to observe, the first mechanism because the frequencies 
of the two lasers, while very stable, may differ by gigahertz. 
If light from one interferometer scatters into the main beam of 
the other, it will likely be at a very different frequency and will 
not produce signals in our 8 kHz band when it beats against 
the reference light for that interferometer. 

Nevertheless, we checked for a correlation produced by 
light from one detector entering the other by looking for the 
calibration signals OJ] injected into one detector in the signal 
of the second detector. During S5, the following calibration 
line frequencies were injected into HI and H2: 46.70 Hz, 
393.10 Hz, 1144.30 Hz (HI) and 54.10 Hz, 407.30 Hz, 
1159.7 Hz (H2). We note here that all those frequencies are 
outside of our analysis bands. We observed no correlation 
beyond the statistical error of the measurement at any of the 
three calibration line frequencies for either of the two detec- 
tors. This check was done for every week and month and for 
the entire S5 data-set. Hence, we conclude that potential sig- 
nals carried by the light in one detector are not coupled into 
the other detector. 

In contrast, we have observed the second scattering mech- 
anism, in which scattered light from the HI beam returns to 
the HI main beam and H2 light returns to the H2 main beam. 
This type of scattering can produce H1-H2 coupling if scat- 
tered light from HI and from H2 both reflect off of the same 


vibrating surface (which modulates the length of the scatter- 
ing paths) before recombining with their original main beams. 
This mechanism is thought to account for the observation that 
shaking the reflective end cap of the 4 km beam tube (just 
beyond an HI end test mass), produced a shaking-frequency 
peak in both HI and H2 GW channels, even though the near- 
est H2 component was 2 km away. However, this scattering 
mechanism is covered by the PEM system since the vibrations 
that modulate the beam path originate in the monitored envi- 
ronment. 

We tested our expectation that scattering-induced correla- 
tions would be identified by our PEM-coherence method. We 
initiated a program to identify the most important scattering 
sites by mounting shakers on the vacuum system at 21 differ- 
ent locations that were selected as potential scattering sites, 
and searching for the shaking signal in the GW channels. All 
significant scattering sites that we found in this way were 
well-monitored by the PEM system. At the site that pro- 
duced the greatest coherence between the two detectors (a re- 
flective flange close to and perpendicular to the beam paths 
of both interferometers), we mounted an accelerometer and 
found that the coherence between this accelerometer and the 
two GW channels was no greater than that for the sensors in 
the pre-existing sensor system. These results suggest that the 
PEM system adequately monitored scattering coupling. As 
we shall show in Sec. IVI Al below. no correlated noise (either 
environmental or instrumental, either narrow-band or broad- 
band) that is not adequately covered by the PEM system is 
identified in the high-frequency analysis, further solidifying 
the adequacy of PEM system. 

V. ANALYSIS PROCEDURE 

In the previous section we described a number of methods 
for identifying correlated noise when searching for a SGWB. 
Here we enumerate the steps for selecting the time segments 
and frequency bands that were subsequently used for the anal- 
ysis. 

STEP 1: We begin by selecting time periods that pass a 
number of data quality flags. In particular, we reject periods 
when: (i) there are problems with the calibration of the data; 

(ii) the interferometers are within 30 s of loss of servo control; 

(iii) there are artificial signals inserted into the data for calibra- 
tion and characterization purposes; (iv) there are PEM noise 
injections; (v) various data acquisition overflows are observed; 
or (vi) there is missing data. With these cuts, the intersection 
of the HI and H2 analyzable time was ~ 462 days for the S5 
run. 

STEP 2: After selecting suitable data segments, we make a 
first pass at determining the frequency bins to use in the anal- 
ysis by calculating the overall coherence between the detector 
outputs as described in Sec. HV A1 Excess coherence levels led 
us to reject the frequencies 86-90 Hz, 100 Hz, 102-126 Hz, 
140.25-141.25 Hz, and 150 Hz in the low-frequency band (80- 
160 Hz), as well as ±2 Hz around the 60 Hz power-line har- 
monics and the violin-mode resonances at 688.5 ± 2.8 Hz 
and 697 ± 3.1 Hz in the high-frequency band (460-1000 
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Hz). It also identified a period of about 17 days in June 
2007 (between GPS times 866526322 and 867670285), dur- 
ing which the detector H2 suffered from excessive transient 
noise glitches. We reject that period from the analysis. 

STEP 3: We perform a search for transient excess power in 
the data using the wavelet-based Kleine Welle algorithm | |45jl , 
which was originally designed for detecting GW bursts. This 
algorithm is applied to the output of both detectors, producing 
a list of triggers for each detector. We then search the two trig- 
ger lists and reject any segment that contains transients with 
Kleine Welle significance larger than 50 in either of the two 
detectors. The value of 50 is a conservative threshold, chosen 
based on other studies done on the distribution of such triggers 
in S5 iH. 

STEP 4: Having determined the reasonably good fre- 
quency bands, we then calculate Cl a and its uncertainty ay, 
summed over the whole band, cf. Eq. [7] The purpose of this 
calculation is to perform another level of data-quality selec- 
tion in the time-domain by identifying noisy segments of 60 s 
duration. It is similar to the non-stationarity cut used in the 
previous analyses 12614281 [46 1 where we remove time seg- 
ments whose a ^ differs, by a pre-determined amount, from 
that calculated by averaging over two neighboring segments. 
Here we use a 20% threshold on the difference. The combi- 
nation of the time-domain data quality cuts described in Steps 
1-4 removed about 22% of the available S5 H1-H2 data. 

STEP 5: After identifying and rejecting noisy time seg- 
ments and frequency bins using Steps 1 — 4, we then use 
the time-shift and the PEM-coherence methods described in 
Secs. IIV Bl and HVCl to identify any remaining contaminated 
frequency bins. To remove bad frequency bins, we split the 
S5 dataset into week-long periods and for each week, we re- 
ject any frequency bin for which either SNRo, y ,ts (/) | or 
|SNRn Q ,pEM(/)| exceeds a pre-determined threshold in the 
given week, the corresponding month, or in the entire S5 
dataset. This procedure generates (different) sets of frequency 
notchings for each week of the S5 dataset. In the analysis 
we use two different sets of SNR threshold values for the cut, 
which are further described in Sec. lVIl 

Figure[4]is a spectrogram of SNRo 0 pem for the 80-160 Hz 
band for all weeks in S5; the visible structure represents cor- 
related noise between HI and H2, which was identified and 
subsequently excluded from the analysis by the H1-H2 coher- 
ence, time-shift, and PEM-coherence measurements. 

Note that previous stochastic analyses using LIGO data 
j26l428[ l46h followed only steps 1, 2 and 4. Steps 3 and 5 
were developed for this particular analysis. 

Having defined the time-segments and frequency-bins to be 
rejected in each week of the S5 data, we proceed with the 
calculation of the estimators and standard errors, Cl a (f) and 
ay, (/), in much the same manner as in previous searches for 
isotropic stochastic backgrounds 112614281 14711 . The data is di- 
vided into T = 60 s segments, decimated to 1024 Hz for the 
low-frequency analysis and 4096 Hz for the high-frequency 
analysis, and high-pass filtered with a 6 th order Butterworth 
filter with 32 Hz knee frequency. Each analysis segment is 
Hann-windowed, and to recover the loss of signal-to-noise 


ratio due to Hann-windowing, segments are 50% overlapped. 
Estimators and standard errors for each segment are evaluated 
with a A / = 0.25 Hz frequency resolution, using the fre- 
quency mask of the week to which the segment belongs. A 
weighted average is performed over all segments and all fre- 
quency bins, with inverse variances, as in Eq. [7] but properly 
accounting for overlapping. 

VI. ANALYSIS RESULTS 

The analysis is separated into two parts corresponding to 
searches for SGWBs with spectral index a = 0 and a = 3 
as described in Sec. [Ill] Since the strain output of an interfer- 
ometer due to GWs is S gw (f) oc /“” 3 (see Eq. [3}. the case 
a = 0 is dominated by low frequencies while a = 3 is inde- 
pendent of frequency. Since for a — 3 there is no preferred 
frequency band, and since previous analyses Phil for stochas- 
tic backgrounds with a = 3 considered only high frequencies, 
we also used only high frequencies for the a = 3 case. Thus, 
the two cases of a = 0 and a = 3 correspond to the analysis 
of the low and high-frequency bands, respectively. In this sec- 
tion, we present the results of the analyses in the two different 
frequency bands as defined in Sec. IIV Dl corresponding to the 
two different values of a. 

To illustrate the effect of the various noise removal methods 
described in the previous two sections, we give the results as 
different stages of cuts are applied to the data (see Table [I]). 
The threshold value used at stage III comes from an initial 
study performed using playground data to understand the ef- 
fectiveness of the PEM-coherence method in finding problem- 
atic frequency bins in the H1-H2 analysis, and hence those 
results are considered as blind analysis results. But a post- 
unblinding study showed that we could lower the SNRpem 
threshold to values as low as 0.5 (for low-frequency) and 1 
(for high-frequency), which are used at stage IV. These post- 
blinding results are used in the final upper-limit calculations. 
For threshold values < 0.5 (low-frequency) or < 1 (high- 
frequency), the PEM-coherence contribution, H a pEM, varies 
randomly as the threshold is changed indicating the statistical 
noise limit of the PEM-coherence method. 


A. High-frequency results 

We performed the high-frequency analysis with spectral in- 
dex a = 3, and reference frequency / re f = 900 Hz. Tables HT1 
and Hill summarize the results after applying several stages of 
noise removal as defined in Table [I] Table [IT] applies to the 
full analysis band, 460-1000 Hz; Table Hill gives the results 
for 5 separate sub-bands. The values of the estimator, fR, the 
PEM-coherence contribution to the estimator, I^pem, and 
the statistical uncertainty, ay, ; , are given for each band and 
each stage of noise removal. Also given is the ratio of the 
standard deviation of the values of the inverse Fourier trans- 
form of ST 3 (/) to the statistical uncertainty oy, , which is a 
measure of excess residual correlated noise. In the absence of 
correlated noise, we expect the distribution of data points in 
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FIG. 4. (Color online) Spectrograms displaying the absolute value of SNRq,pbm(/) for 80-160 Hz (left) and 460-1000 Hz (right) as a 
function of the week in S5. The horizontal dark (blue) bands correspond to initial frequency notches as described in STEP 2 (Sec. 0 and 
vertical dark (blue) lines correspond to unavailability of data due to detector downtime. The large SNR structures seen in the plots were 
removed from the low- and high-frequency analyses. 


Stage 

High-frequency analysis 

Low-frequency analysis 

Steps 

% of data 
vetoed 

Steps 

% of data 
vetoed 

I 

Step 1 

8.51 

Step 1 

8.51 

II 

Steps 1-4 

35.88 

Steps 1^1 

56.01 

III 

Steps 1-5 with 
|SNRpem| > 2 , 
|SNRts| > 2 

47.19 

Steps 1-5 with 
|SNRpem| > 2 , 
|SNRts| > 2 

72.29 

IV 

Steps 1-5, with 
|SNRpem| > 1, 
|SNRts| > 2 

48.95 

Steps 1-5, with 
|SNRpem| > 0.5, 
|SNRts| > 2 

76.60 


TABLE I. Definition of various stages of noise removal for the high 
and low-frequency analyses in terms of the analysis steps described 
in Sec. 0 Here stage III corresponds to the blind analysis and stage 
IV to the post-unblinding analysis. The percentage of data vetoed 
accounts for both the time segments and frequency bins excluded 
from the analysis. In calculating veto percentage, the analyses with 
non-colocated LIGO detectors only accounts for the time segments 
excluded from the analyses and is the reason for the large numbers 
we see in the last column compared to other LIGO analyses. 


the inverse Fourier transform of O 3 (/) to follow a Gaussian 
distribution with mean 0 and std . Hence a ratio 1 is 
a sign of excess correlated noise, which shows up as visible 
structure in the plot of the inverse Fourier transform of ( 13 (f) 
(for example, see the right hand plots in Fig. [3. We see that 
this ratio decreases for the full 460-1000 Hz band and for each 
sub-band with every stage of data cleanup indicating the ef- 
fectiveness of PEM-coherence SNR cut. We also note that the 
values listed in Tables UTI mil and HVl are the zero lag values of 
Cl n in the corresponding inverse Fourier transform plots. 

Figure 0 is devoted entirely to the noisiest sub-band, 


Stage 

CI 3 

(xlO -4 ) 

H3.PEM 

(xlO -4 ) 

°n 3 

(xlO -4 ) 

std/o-fi 3 

I 

77.5 

—3.05+ 

2.82 

20.5 

II 

-2.17 

-3.62 

3.24 

1.18 

III 

-4.11 

-4.30 

3.59 

1.04 

IV 

-1.29 

-2.38 

3.64 

1.01 


TABLE II. Results for the H1-H2 high-frequency analysis (460- 
1000 Hz) after various stages of noise removal were applied to the 
data. The estimates Cl 3 , PEM-coherence contribution, H 3 ,pem and 
<Tq 3 are calculated assuming Ho = 68 km/s/Mpc. < 7 ^ is the sta- 
tistical uncertainty in CI 3 . The last column gives the ratio of the stan- 
dard deviation of the values of the inverse Fourier transform of CI 3 (/) 
to the statistical uncertainty Cq 3 - As described in Sec. IVI Al a ratio 
much 1 is a sign of excess cross-correlated noise. Mhe PEM- 
coherence estimate on stage I also excludes frequencies (including 
60 Hz harmonics) and time segments similar to stages II-IV. 


628-733 Hz. The left column of plots shows S^ 3 (/) and 
^ 3 .pem(/)i with black lines denoting the statistical error bar 
±tj^ 3 (/). Here we can clearly see the effectiveness of noise 
removal through the four stages discussed above. Note the 
lack of structure near zero-lag in the final inverse Fourier 
transform of the estimator S 2 3 (/) which is consistent with no 
correlated noise. Figure [ 6 ] is a similar plot for the full 460- 
1000 Hz band, showing the results after the final stage of cuts. 
Again note the lack of significant structure near zero-lag in 
the inverse Fourier transform of ( 13 (f). Figure [7] (left panel) 
shows how the final estimate, CI 3 , summed over the whole 
band, evolves over the course of the run after the final stage 
of cuts. The smoothness of that plot (absence of any sharp 
rise or fall after the accumulation of sufficient data i.e., one 
month) indicates that no particular time period dominates our 
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Band 

Stage H3 

H3.PEM 

°A 3 std/crn 3 

(Hz) 


(xlO -4 

) (xlO -4 ) (xlO -4 ) 


460-537 

I 

-7.28 

-0.22 

4.48 

5.40 


II 

-2.17 

-0.24 

5.08 

1.01 


III 

-0.60 

-1.23 

5.68 

0.98 


IV 

-0.34 

-1.23 

5.69 

0.97 

537-628 

I 

163 

-2.28 

5.46 

24.0 


II 

14.7 

-2.46 

6.32 

1.08 


III 

8.83 

-2.00 

6.96 

1.02 


IV 

8.56 

-1.98 

7.03 

1.02 

628-733 

I 

512 

-16.7 

7.33 

35.9 


II 

-33.2 

-20.5 

8.52 

1.37 


III 

-37.0 

-16.3 

9.20 

1.21 


IV 

-26.5 

-5.88 

9.66 

1.12 

733-856 

I 

-397 

-1.77 

8.32 

23.0 


II 

-4.44 

-2.24 

9.49 

1.67 


III 

-5.29 

-6.40 

11.0 

1.04 


IV 

2.76 

-3.91 

11.3 

0.98 

856-1000 

1 

89.2 

4.63 

10.6 

3.37 


II 

2.44 

4.63 

12.0 

1.02 


III 

0.004 

-1.47 

13.2 

1.01 


IV 

0.21 

-1.41 

13.2 

1.01 

TABLE III. Same as 

Table ITT1 but for 5 separate sub-bands of 

1000 Hz. 






Stage 


Ho fio.pEM CT n 0 

std/crn. 

0 


(xlO -6 ) (xlO -6 ) (xl0“ 

6 ) 


I 

6.17 —0.39* 0.44 

5.90 


II 

- 

1.71 - 

0.78 0.63 

1.80 


III 

- 

1.57 - 

0.84 0.79 

1.64 


IV 

- 

0.26 - 

0.29 0.85 

1.63 



TABLE IV. Similar to Tablelillbut for the low-frequency analysis (80— 
160 Hz) and for spectral index a = 0. The different rows give the 
results after various stages of noise removal were applied to the data. 
^The PEM-coherence estimate on stage I also excludes frequencies 
(including 60 Hz harmonics) and time segments similar to stages II- 
IV. 


final result. 


B. Low frequency results 

We now repeat the analysis of the previous subsection but 
for the low-frequency band, 80-160 Hz with spectral index 
a = 0 and / re f = 100 Hz. Table HVl summarizes the results 
for the low-frequency analysis after applying several stages 
of noise removal as defined in Table U Figure [8] shows the 
results obtained by applying the noise removal cuts in four 
stages. The left column of plots contain the estimators, fig(/) 
and Ho,pem(/), with lines denoting the statistical error bar 
±fJ n 0 (/)- 


In contrast to the high-frequency analysis (compare Figs. [6] 
and [8} there is still much structure in the inverse Fourier 
transform of Og(f) around zero-lag even after the final stage 
of noise removal cuts were applied. In addition, the PEM- 
coherence contribution to the estimator, Ho.pem (/), displays 
much of the structure observed in 0q(/). Both of these ob- 
servations suggest contamination from residual correlated in- 
strumental or environmental noise that was not excluded by 
the noise removal methods. Figure [7] (right panel) shows how 
the final estimate, T2o, evolves over the course of the run after 
the final stage of cuts. We note here that even though Ho (last 
entry in Table II Vt is consistent with zero (within 2a), its esti- 
mate at other non-zero lags vary strongly as shown in Fig. [8] 
(lower right). This indicates the presence of residual corre- 
lated noise after all the time-shift and PEM-coherence noise 
removal cuts are applied. 


C. Hardware and software injections 

We validated our analysis procedure by injecting simulated 
stochastic GW signals into the strain data of the two detec- 
tors. Both hardware and software injections were performed. 
Hardware injections are performed by physically moving the 
interferometer mirrors coherently between interferometers. In 
this case the artificial signals were limited to short durations 
and relatively large amplitudes. The data from these hardware 
injection times were excluded from the analyses described 
above, as noted in Sec.[V] Step 1. Software injections are con- 
ducted by adding a simulated GW signal to the interferometer 
data, in which case they could be long in duration and rela- 
tively weak in amplitude. During S5 there was one stochastic 
signal hardware injection when both HI and H2 were oper- 
ating in coincidence. A stochastic background signal with 
spectral index a = 0 and amplitude Ho = 6.56 x 10~ 3 
was injected for approximately 3 hours. In performing the 
analysis, frequency bins were excluded based on the standard 
H1-H2 coherence calculations. No additional frequency bins 
were removed using SNRpem- The recovered signal was 
Ho = (7.39 ± 1.1) x 10~ 3 , which is consistent with the 
injected amplitude. Due to the spectral index used for the 
injection (a = 0), the recovery analysis was performed us- 
ing only the low frequency band. We also performed a soft- 
ware injection in the high frequency band with an amplitude 
H3 = 5.6 x 10 -3 , and we recovered it successfully. Figure[9] 
shows the spectrum of the recovered 0 . 3 (f) and its inverse 
Fourier transform. 


VII. ASSESSING THE RESIDUAL CORRELATED NOISE 

After applying the full noise removal procedure, the high- 
frequency band appears clean whereas the low-frequency 
band exhibits evidence of residual correlated noise. In order 
to interpret the implications of these two very different results, 
we introduce a general procedure for determining whether 
a stochastic measurement is sufficiently well-understood to 
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FIG. 5. Plots of 03 (f) and H 3 ,pem (/) (left), and the inverse Fourier transform of 0 . 3 (f) (right), for the (noisiest) 628-733 Hz sub-band after 
various stages of noise removal were applied to the data. The four rows correspond to the four different stages of cleaning defined in Tabic [T] 
(The top right plot has y-axis limits 13 x greater than the other three.) 
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FIG. 6. Plots of ( 13 (f) and fi 3 ,pEM(/) (left), and the inverse Fourier transform of complex Cls(f) (right) for the full band (460-1000 Hz) 
after the final stage of noise removal cuts. 



FIG. 7. Running point estimates CI 3 and Clo for the high-frequency (460-1000 Hz) and low-frequency (80-160 Hz) analyses, respectively (left 
and right panels). The final stage of noise removal cuts have been applied for both analyses. 


yield an astrophysical interpretation. While our immediate 
concern is to provide a framework for interpreting the two 
results presented here, we aim to give a comprehensive pro- 
cedure that can be applied generally, to both co-located and 
non-co-located detectors. In this spirit, this section is orga- 
nized as follows: first, we present a general framework for 
interpreting stochastic measurements; then we discuss how it 
can be applied to (familiar) results from non-co-located detec- 
tors; and finally we apply the framework to our present results. 

To determine whether a result can be interpreted as a con- 
straint on the SGWB, we consider the following three criteria: 

1 . We have accounted for all known noise sources through 
either direct subtraction, vetoing, and/or proper estima- 
tion of systematic errors. 

2. Having accounted for known noise sources, we do not 
observe evidence of residual noise that is inconsistent 
with our signal and noise models. 


3. To the best of our knowledge, there is no plausible 
mechanism by which broadband correlated noise might 
be lurking beneath the uncorrelated noise at a level com- 
parable to the GW signal we are trying to measure. 

If an analysis result does not meet these criteria, then we con- 
servatively place a bound on the sum of the GW signal and the 
residual correlated noise. If a result meets all the criteria, then 
we present astrophysical bounds on just the GW signal. 

Let us now examine these criteria in the context of previous 
results using the non-co-located LHO and LLO detectors |j28tl . 
Criterion #1 was satisfied by identifying and removing instru- 
mental lines attributable to known instrumental artifacts such 
as power lines and violin resonances. Criterion #2 was satis- 
fied by creating diagnostic plots, e.g., showing S2o vs. lag (the 
delay time between the detectors; see Fig. 0. which demon- 
strated that the measurement was consistent with uncorrelated 
noise (and no GW signal). Criterion #3 was satisfied by per- 
forming order-of-magnitude calculations for plausible sources 
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FIG. 8. Plots of Clo(f) and Ho.pem (/) (left), and the inverse Fourier transform of Qo (/) (right) for the 80-160 Hz band after various stages 
of noise removal were applied to the data. The four rows correspond to the four different stages of cleaning defined in Tabled 
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FIG. 9. Left panel: Recovered spectrum for a software injection with an amplitude fl 3 = 5.6 x 10 -3 (SNR ~ 17). Right panel: The inverse 
Fourier transform of the recovered Cis (/) and a ±10ms zoom-in around zero lag. 


of correlated noise for LHO-LLO including electromagnetic 
phenomena, and finding that they were too small to create 
broadband correlated noise at a level that is important for ini- 
tial LIGO. 

Next, we consider how the criteria might be applied to fu- 
ture measurements with non-co-located detectors. During the 
advanced detector era, correlated noise from Schumann reso- 
nances may constitute a source of correlated noise at low fre- 
quencies < 200 Hz, even for widely separated detector pairs 
such as LHO-LLO ll29L [30h . While it may be possible to mit- 
igate this potential correlated noise source through commis- 
sioning of the detectors to minimize magnetic coupling, or 
failing that, through a noise subtraction scheme, we consider 
the possibility that residual correlated noise is observed. In 
this scenario, we could still aim to satisfy criteria #2 and #3 by 
using magnetometer measurements to construct a correlated 
noise budget, which could then be used to interpret the results. 

Finally, we consider how the criteria apply to the measure- 
ments presented in this paper. The high-frequency analysis 
meets criteria 1 and 2 as we did not observe residual noise in- 
consistent with our noise models (see Fig. [6}. We did observe 
residual noise for the low frequency analysis (see Fig. [8}, but 
it was consistent with a preliminary noise model, based on 
measured acoustic coupling and microphone signals (most of 
the channels identified by the PEM coherence method were 
either microphones or accelerometers placed on optical tables 
that were susceptible to acoustic couplings). While the bands 
that were acoustically loudest (containing certain electronics 
fans) were vetoed, the acoustic coupling in between the ve- 
toed bands was high enough to produce a residual signal. We 
did not further develop the noise model to meet criterion 1 be- 
cause, with the systematic error from acoustic coupling, the 
astrophysical limit would not have improved on values we 
have reported previously [28, 48]. For this reason, we do not 
present an astrophysical limit for the low frequency band. 

We addressed criterion #3 in two ways. First, by investigat- 
ing mechanisms that might produce un-monitored broad-band 
correlations between detectors, such as the study of correla- 


tions introduced by the shared data acquisition system, the 
study of correlations introduced by light scattering, and PEM 
coverage studies described in Sec. lIV El 

We also identified the sources of most of the features be- 
tween 80 and 400 Hz. For many of the spectral peaks, in 
addition to coherence between the GW channels, there was 
also coherence between the individual GW channels and the 
accelerometer and microphone signals from the vertex area 
shared by both detectors. The coupling was consistent with 
the measured coupling of acoustic signals to the detectors. 
Most of these features were traced to electronics cooling fans 
in specific power supply racks in the vertex station by compar- 
ing coherence spectra to spectra for accelerometers mounted 
temporarily on each of the electronics racks. The features 
were produced at harmonics of the fan rotation frequencies. 

The second type of coherence feature was associated with 
bilinear coupling of low frequency (< 15 Hz) seismic motion 
and harmonics of 60 Hz, producing side-band features around 
the harmonics that were similar to the features in the 0-15 Hz 
seismic band. Coherence of side-band features was expected 
since the coherence length of low-frequency seismic signals 
was greater than the distance separating sensitive parts of the 
two interferometers at the vertex station, and the seismic iso- 
lation of the interferometers was minimal below 10 Hz. 

In conclusion, we found no peaks or features in the coher- 
ence spectrum for the two GW channels that were inconsis- 
tent with linear acoustic coupling or bilinear coupling of low 
frequency seismic noise and 60 Hz harmonics at the vertex 
station. Neither of these mechanisms is capable of producing 
broad-band coherence that is not well monitored by the PEM 
system. Therefore, for the high frequency analysis, we satisfy 
the three criteria for presenting astrophysical bounds on just 
the GW signal. 
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A. Upper-limits 

Since there is no evidence of significant residual noise con- 
taminating the high-frequency data after applying the full 
set of cuts, we set a 95% confidence-level Bayesian upper- 
limit on fl3. We use the previous high-frequency upper limit 
f^3 < 0.35 (adjusted for // 0 = 0.68 km/s/Mpc) from the 
LIGO S5 and Virgo VSR1 analysis 14(^ 1 as a prior and assume 
a flat distribution for 0 : > from 0 to 0.35. We also marginal- 
ize over the calibration uncertainty for the individual detectors 
(10.2% and 10.3% forHl and H2, respectively). In order to in- 
clude in our calculation the PEM estimate of residual contam- 
ination, we take cr? + fig PEM as our total variance. We note 

here that the estimated is within the observed a {) 

i.e., we observe no evidence of excess environmental contami- 
nation and the above quadrature addition increases the limit by 
~ 20%. The final result is Qg < 7.7 x 10 -4 for the frequency 
band 460-1000 Hz, which is an improvement by a factor of 
~ 180 over the recent S6/VSR2-3 result ll48ll . All of the above 
~ 180 factor improvement comes from the nearly-unity over- 
lap reduction function of the co-located Hanford detectors. In 
fact, all other data being same, if we were to consider the H2 
detector to not be located at Hanford but instead at the LIGO 
Livingston site yields an upper limit that is worse by a factor 
of ~ 1.7 than the S6/VSR2-3 result. Most of this difference 
of ~ 1.7 comes from the improved sensitivities of S6/VSR2-3 
detectors compared to S5 H1-H2 detectors. Upper-limits for 
the five separate sub-bands of the high-frequency analysis are 
given in Table [V] 


Band (Hz) 

95% CL UL (xlO -3 ) 

460-1000 

0.77 

460-537 

1.11 

537-628 

2.12 

628-733 

1.18 

733-856 

2.53 

856-1000 

2.61 


TABLE V. 95% confidence level upper-limits for the the full band 
(460-1000 Hz) and for five separate sub-bands. 

As mentioned in Sec. IVIBI the structure in the inverse 
Fourier transform plots of Fig.[8]suggests contamination from 
residual correlated noise for the low-frequency analysis and 
hence we do not set any upper-limit on Ho using the low- 
frequency band 80-160 Hz. 

VIII. SUMMARY AND PLANS FOR FUTURE ANALYSES 

In this paper, we described an analysis for a SGWB using 
data taken by the two co-located LIGO Hanford detectors, HI 
and H2, during LIGO’s fifth science run. Since these detectors 
share the same local environment, it was necessary to account 
for the presence of correlated instrumental and environmental 
noise. We applied several noise identification and mitigation 
techniques to reduce contamination and to estimate the bias 


due to any residual correlated noise. The methods proved to 
be useful in cleaning the high-frequency band, but not enough 
in the low-frequency band. 

In the 80 — 160 Hz band, we were unable to sufficiently 
mitigate the effects of correlated noise, and hence we did not 
set any limit on the GW energy density for a = 0. For the 
460 — 1000 Hz band, we were able mitigate the effects of cor- 
related noise, and so we placed a 95% C.L. upper limit on the 
GW energy density alone in this band of H3 < 7.7 x 10 -4 . 
This limit improves on the previous best limit in the high- 
frequency band by a factor of ~ 180 148]] . Figure flOl shows 
upper limits from current/past SGWB analyses, as well as lim- 
its from various SGWB models, and projected limits using 
Advanced LIGO. We note here that the indirect limits from 
BBN apply to SGWBs present in the early universe at the 
time of BBN (and characterized by an a = 0 power law; see 
Eq.0, but not to SGWBs of astrophysical origin created more 
recently (and assumed to be characterized by an a = 3 power 
law). Thus, the results presented here complement the indi- 
rect bound from BBN, which is only sensitive to cosmologi- 
cal SGWBs from the early universe, as well as direct a = 0 
measurements using lower-frequency observation bands [28(1 . 

There are several ways in which the methods presented in 
this paper can be improved. We list some ideas below: 

(i) As mentioned in Sec. lIV Cl we can improve the estimate 
of the PEM contribution to the coherence by allowing for cor- 
relations between different PEM channels zi and zj. This 
requires inverting the full matrix of PEM coherences 7 jj(/) 
or solving a large number of simultaneous equations involving 
7zj(/), rather than simply taking the maximum of the prod- 
uct of the coherences as was done here. A computationally- 
cheaper alternative might be to invert a sub-matrix formed 
from the largest PEM contributors — i.e., those PEM channels 
that contribute the most to the coherence. 

(ii) We can use bicoherence techniques to account for (non- 
linear) up-conversion processes missed by standard coherence 
calculations. This may allow us to identify cases where low- 
frequency disturbances excite higher-frequency modes in the 
detector. 

(iii) The estimators Cl a (f) used in this analysis are optimal 

in the absence of correlated noise. In the presence of corre- 
lated noise, these estimators are biased, with expected values 
given by the sum, + r] a (f), where the second term in- 
volves the cross-spectrum, of the noise contribution 

to the detector output. An alternative approach is to start with 
a likelihood function for the detector output si, § 2 , where we 
allow (at the outset) for the presence of cross-correlated noise. 
(This would show up in the covariance matrix for a multi- 
variate Gaussian distribution.) We can parametrize Niz(f) in 
terms of its amplitude, spectral index, etc., and then construct 
posterior distributions for these parameters along with the am- 
plitude and spectral index of the stochastic GW signal. In this 
(Bayesian) approach, the cross-correlated noise is treated on 
the same footing as the stochastic GW and is estimated (via its 
posterior distribution) as part of the analysis I5~8h . However, 
as described in |5j|, this works only for those cases where the 
spectral shapes of the noise and signal are different from one 
another. 
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FIG. 10. (Color online) Upper limits from the current H1-H2 analy- 
sis, previous SGWB analyses and the projected advanced LIGO limit, 
along with various SGWB models.The BBN limit is an integral limit 
on fl gw i.e., f in the 10~ 1 2 3 4 5 6 7 8 9 10 — 10 10 Hz band derived 

from the Big Bang Nucleosynthesis and observations of the abun- 
dance of light nuclei f28l,[491. The measurements of CMB and matter 
power spectra provide a similar integral bound in the frequency range 
of 10~ 15 — 10 10 Hz(H. The pulsar limit is a bound on the fl gw (/) 
at / = 2.8 nHZ and is based on the fluctuations in the pulse arrival 
times from millisecond pulsars m. In the above figure, slow-roll 
inflationary model f52fl assumes a tensor-to-scalar ratio of r = 0.2, 
the best fit value front the BICEP2 analysis d- In the axion based 
inflationary model, for certain ranges of parameters the backreaction 
during the final stages of inflation is expected to produce strong GWs 
at high frequencies |[53j]. The stiff equation of state (EOS) limit cor- 
responds to scenarios in the early universe (prior to BBN) in which 
GWs are produced by an unknown ‘stiff’ energy |54tl . For the above 
figure we used the equation of state parameter w = 0.6 in stiff EOS 
model. The cosmic string model corresponds to GWs produced by 
cosmic strings in the early universe (55}] . The Earth’s normal mode 
limits are based on the observed fluctuations in the amplitudes of 
Earth’s normal modes using an array of seismometers [56]]. The as- 
trophysical SGWBs (BBH and BNS) are due to the superposition of 
coalescence GW signals from a large number of binary black holes 
(BBH) and binary neutron stars (BNS) (U. 


(iv) We can also reduce correlated noise by first removing 
as much noise as possible from the output of the individual de- 
tectors. Wiener filtering techniques can be applied to remove 
acoustic, magnetic, and gravity-gradient noise from the time- 
series output of the LIGO detectors ll60l - l62ll . Furthermore, 
feed-forward control can be used to to cancel seismically- 
induced motion before it affects the LIGO test masses M- 


These and/or other techniques might be needed for future 
cross-correlation searches using advanced detectors, where 
improved (single-detector) sensitivity will mean that corre- 
lated noise may be an issue even for physically-separated de- 
tectors, such as the LIGO Hanford-LIGO Livingston detector 
pair 
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